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ABSTRACT 



> 

Q ■ Summation over hard thermal loops, by themselves and as insertions in higher 

order Feynman diagrams, is important in thermal perturbation theory for Quantum 

Chromodynamics, so that all contributions of a given order in the coupling constant 

5v ', can be consistently taken into account. I review some of the basic properties of hard 

thermal loops and how the generating functional for them is related to the eikonal 

for a Chern-Simons gauge theory, and using an auxiliary field, to the gauged Wess- 

^ ■ Zumino-Novikov-Witten action. The Hamiltonian analysis of the effective action 

^D \ and a discussion of plasma waves are given. It is also pointed out that a possible 

^^ ■ expression for the magnetic mass term can be written in a closely related way. 
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1. Introduction 

In these lectures, I review ttie relationstiip between Cliern-Sinions (CS) and 
Wess-Zumino-Novikov-Witten (WZNW) ttieories and Quantum Ctironiodynaniics 
(QCD) at finite temperature ^~^. Tiie results are expected to be useful in under- 
standing tlie tlieory of the quark-gluon plasma. 

There is essentially universal agreement by now that QCD is the theory of 
strong interactions. One of the interesting features of QCD is the deconfinement 
phase transition which is indicated by general theoretical arguments as well as 
lattice simulations ^'^. If the ambient temperature is raised beyond a critical value 
Tc, hadronic matter makes a phase transition. Quarks and gluons are no longer 
confined; we have a plasma of quarks, antiquarks and gluons. As a new phase of 
hadronic matter, the quark-gluon plasma is of considerable interest to physicists. 

The plasma phase of hadronic matter, it is generally believed, can be achieved 
in the collisions of sufficiently heavy nuclei at sufficiently high energy. However, the 
key word here is sufficiently. Will the relativistic heavy ion collider (RHIC), to be 
built at Brookhaven, have enough energy to achieve this? Can thermal equilibrium 
be obtained even for a short duration before the plasma cools and rehadronizes? 
These are questions to which, as of now, there are no definitive and clear answers ^. 
Nevertheless, understanding physical phenomena in the quark-gluon plasma at and 
near thermal equilibrium, which can be described by QCD at high temperatures, is 
an important and necessary first step. The critical temperature Tc for the transition 
to the plasma phase is expected to be of the order of the scale parameter, Aqcd 
which is around 200 MeV. Physical phenomena near this temperature are of great 
interest, especially since they may serve as signatures for the transition '^. We 
shall not discuss these however. Instead, we shall consider temperatures which are 
significantly higher than Tc] the average momentum carried by quarks and gluons 
will be much higher than Aqcd and one can expect thermal perturbation theory to 
provide an adequate description of the physics. We are interested in the technical 
problems of thermal perturbation theory for QCD, especially hard thermal loops and 
related physical phenomena such as Debye screening. Landau damping, propagation 
of plasma waves, etc. (Eventhough we consider T > Tc, many of the results are of 
importance near Tc as well. It is qualitatively easier to isolate the relevant terms 
at high temperatures; the modifications required at lower temperatures can be 
easily included. The quark flavors of interest are also the light quarks, with masses 
small compared to T. The heavier flavors have the standard Boltzmann suppression 
factors and lead to small modiflcations.) 

As is well known , among the most elementary phenomena in a plasma are Debye 
screening and Landau damping. If we consider an Abelian plasma of positive and 
negative charges ±e, viz., electrodynamics, screening can be understood using the 
classic argument of Debye. One considers the Poisson equation for the electrostatic 
potential of a test charge, say positive, in the plasma. 

'geAo/T _ g-eAo/T^ 
eAo/T 1 p-eAo/T 
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where the right hand side is the charge density in the vicinity of the test charge. 
n is the average number density of particles. The exponentials are the Boltzmann 
factors giving the preferential accumulation of negative and depletion of positive 
charges in the vicinity due to the Coulomb forces, with proper normalization. The 
approximation (lb), which is valid for high temperatures, shows that the solutions 
have the screened Coulomb form (l/r) exp(-TOc?') with a Debye screening mass m|, = 
(ne^/T). For a relativistic plasma, the qualitative features of this argument are valid 
and with n ^ T^, we expect to|, ^ e^T^. And by calculating the photon propagator 
in thermal electrodynamics, one can indeed obtain a similar result *. 

The above argument is simple and nice, but is presented in terms of potentials 
and as such, it does not seem to be gauge invariant. For the Abelian case, one 
can easily reformulate the arguments using only the gauge-invariant electric and 
magnetic fields. However, in a non- Abelian plasma, such as the quark-gluon plasma, 
it is difficult to avoid the use of gauge potentials altogether. Further, since even 
the notion of the charge of a gluon has to be defined with respect to some chosen 
Abelian direction of the gauge group, it is clear that the simple argument of Eq.(l) 
will have to be modified. The question of interest to us is: how do we obtain a 
gauge-invariant description of Debye screening in a non- Abelian plasma (in terms 
of gauge potentials)? More specifically, we need a functional of the gauge potentials, 
r[yl], which is generated by the statistical distributions and is effectively a gauge- 
invariant mass term for the gauge fields. Screening, of course, is the static part of 
a more general problem, the dynamical part of which is the propagation of plasma 
waves. Such a T[A] will thus be important for the discussion of plasma waves as 
well. 

Landau damping arises from the scattering of the particles in the plasma against 
an external field. This can transfer energy and momentum from the field to the 
particle, leading to damping of the field. The process is essentially the inverse of 
Ccrenkov radiation and can be described by the proper continuation of the results for 
plasma waves to a region of spacelike momenta. Of course, being a time-asymmetric 
process, one must describe it in terms of the equations of motion rather than an 
action; the relevant current for these equations will be a continuation of what is 
given by T[A]. 

Hard thermal loops are closely related to the above. They arise because of 
the need to carry out a partial resummation of perturbation theory in thermal 
QCD. The need for resummation is easily seen using the following argument due to 
Pisarski ^. Consider the elementary polarization diagram of gluons, fig. (la) and its 
first correction fig. (lb). Before the final loop integration over p, the ratio of diagram 
(lb) to (la) is n(p)/p2. Because of Debye screening, the polarization tensor n(p) has 
a small p-expansion, n(p) ^ g'^T^ + g'^T\p\ + .... {g is the quark-gluon or gluon-gluon 
coupling constant.) We thus see that for the small p-regime of integration, i.e., 
p^ < g^T^, the naively higher order diagram (lb) is comparable to the lower order 
term (la). Therefore, to be consistent to a given order in g, one must sum over 



(la), (lb), and a series of further insertions of n(p). This resummation leads to new 
effective propagators; in general, we shall need new effective vertices as well. The 
resummation can be summarized by saying that the new propagators and vertices 
arise from an action 

S ^ [d^x {-\F') + r[A] (2) 



where we add a term r[^] to the standard Yang-Mills action for gluons. 

Now, going back to figures (1), we see that the regime of interest is when the 
momentum external to n(p), viz., p is of the order gT or less. We can thus take the 
momenta carried by the potentials in r[^] to be ^gT. In order to determine which 
diagrams can contribute to T[A], we shall need to extend the argument given above 
for n(p) to general vertices and obtain rules for counting powers of g and T, keeping 
in mind the restriction on external momenta. These rules basically follow from the 
thermal propagator and let us briefly examine them. 

2. Analysis of diagrams and power counting rules 

The propagator for a quark at finite temperature is defined by 

S{x,y) = {T qix)qiy)) (3) 

where T denotes time-ordering. We shall use Minkowski space propagators with 
thermal averages for the products of creation and annihilation operators; this is 
conceptually the simplest for us since we need to talk about nonequilibrium phe- 
nomena such as the propagation of plasma waves. The propagator can be evaluated 
as 

S{x,y) = I (03 ^{e(^° - y')[a,e-^P-(^-y^-/ ■? + 0,e^P'<^-yh -P']- 
Q(y" - x°)[(3pe-'P<'=-y'>j ■ p + ape'P'<''-y^j ■ p']} 
where p° = \p\, p= {p^,p), p' = {p^,-p), and 6(a;), of course, is the step function. Also 

Op = 1 - np, I3p = Hp. (5) 

The distribution functions np,np are defined by the thermal averages 

(at"''-a^^") = np(5''^(5"'3 

(4"'''^p'') = np5''''5"^ (6) 

where (aj"''^,ap'''), (5j"''',6p''') are the annihilation and creation operators for quarks 
and antiquarks respectively, a, (3 are spin indices; r, s are color indices. For a plasma 
of zero fermion number, we can take 



For a plasma with a nonzero value of fermion number, there is a chemical potential 



and correspondingly 
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are not equal. 



We can also write the propagator, for zero chemical potential, as 



S{x,y) = 
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The first term in brackets is the standard Feynman term at T = 0. The T-dependent 
correction is effectively on mass-shell, enforced by the (5-function. 

For the gluon field, the propagator in the Feynman gauge can be similarly written 
as G'^'l{x,y) = d'^''g^„G{x,y) where g^^ is the metric tensor, a,b are color indices and 



G{x,y) = 



(Pp 1 



[e(x° - y") {ape-'P(^-^) +/3pe'P(^-^)} 



d^p 



(27r) 
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Here ap = 1 + Up, Pp = Up and Up is the bosonic distribution function 

1 

"^P ~ eP"/T - 1 



(9a) 
(96) 

(10) 



For our power counting rules, it is easiest to consider the propagators with the 
time-orderings separated as in Eqs.(4,9a). Suppose now that these propagators are 
in a loop diagram, with integration over momenta. The distributions tell us that 
the average momentum is of the order of T. For each diagram, we thus get ^: 

1) T^ corresponding to d^p 

2) 1/T for each propagator from the l/2p" factor 

3) a power of T for each factor of p^ in the numerators in vertices 

4) for each propagator, other than the first, a power of 1/fe, where k is the external 
momentum (^ gT) 

5) a factor fc/T, one for each loop, for two or more propagators of the same 
statistical type. 

The last two rules can be understood as follows. In using Eqs.(4,9a) for the 
propagators in any Feynman diagram, we encounter many terms corresponding 
to different ordering of the time arguments x°,y°,z°,etc. An obvious strategy for 
simplification is to carry out the time-integrations first, introducing convergence 
factors e"^""^" , e'^'^y" etc., e small and positive, as required. The integrations give 
energy-denominators and bring the result to a form where simplification due to 
the external momenta being soft (^ gT) relative to the loop momenta (^ t) can be 
implemented easily. Among energy-denominators, there will be some of the form 
p° - g" ± A:°, where p,q refer to loop momenta. Because their difference is involved, 
such a denominator is of the order of k, eventhough p°, q" can be of the order 



of T. This gives the 1/fc-factor of rule 4. If we have only one propagator, time- 
integration gives only an energy-conservation (5-function; thus the l/fc-factors are 
only for propagators other than the first. The last rule has to do with the fact 
that the difference of statistical distributions such as Up - Uq is what arises. With 
p- q Ki k, this gives rule 5. The last two rules will become clearer in the course of 
the evaluation of the two-point function which we shall do shortly. 

We can apply these rules, as an example, to the three-point function and its first 
correction as shown. Because of the derivative coupling, fig. (2a) goes as gk ~ g'^T. 
For fig. (2b), first of all we have g^, then T^ from rule 1, l/T^ from rule 2, T^ from 
rule 3, 1/k'^ from rule 4 and k/T from rule 5. With k ^ gT, this is again ^ g'^T. Thus 
fig. (2b) must be included in our calculation of r[A]. Continuing in a similar way, one 
can see that all diagrams which should be included in r are one- loop diagrams ^; 
further, only the T-dependent terms of these one-loop diagrams (and of course with 
external momenta small compared to t) are important. Thermal one-loop diagrams 
with loop momenta ^ T and external momenta ^gT are called hard thermal loops. 

3. Calculation of the two-point function 

Let us now carry out an explicit calculation, say for the two-point function. 
This may be a trifie tedious, but is well worth the effort since it illustrates many of 
the features which generalize to n-point functions ^ We consider a one-loop quark 
graph with two external gluon lines, in other words, the gluon polarization diagram. 
The relevant part of the Lagrangian for the quark fields g, q is 

C^qi-1-{d + A) q (11) 

where A^ = -it'^A'^ is the Lie-algebra-valued gluon vector potential, t°- are hermitian 
matrices corresponding to the generators of the Lie algebra in the representation 
to which the quarks belong; [t°,i''] = if°''"'t'', where /"'"' are the structure constants, 
and Tr(i"t^) = i(5"''. We take the gauge group to be SU{N) with Np fiavors of quarks. 
Prom now on, we shall not explicitly display the quark-gluon coupling constant g, 
as it is easily recovered at any stage hj A —> gA. The one-loop quark graphs are 
given by the effective action 

r = -i Tr log (1 + S* 7 • A) (12) 

The two-gluon term in r is given by 

r(2) = 1 1 d^a:d^y Tr [ 7 • A{x)S{x, y)j ■ A{y)S{y, x)] (13) 

In using Eq.(4) for the propagator, we find four terms in r*^) with a;° > y° and four 
terms with 2/° > x". Writing 

M^) ^ I ^ e""-" Mk) (14) 



and carrying out the time-integrations we get 
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T{p',q'){ 



p'^ + q^-k^-ie p° + q'^ - k'^ + ie' 

apCtq fipPq ^ 

p° + g° + fcO-ie p" + gO + A:0 + ie'^ 

OLpPq fjpCXq 



p"-q° + k^-ie p^ - q" + k^ + ie' 



(15) 



where 



T{p, q) = Tr[ 7 • A(fc)7 • p-y ■ A{k')j ■ q 

d^k d^k' 



dp{k)^{2TTf5^'^\k + k') 



(27r)4 {2ttY 



(16) 

(17) 



In Eq.(15), p = g + k. Since p° = \q + fc| ~ (j° + g • fc/g° for \k\ small compared to |g|, 
the denominators in Eq.(15) involve fc • Q, k- Q' and 2g" + A; • Q, 2g° + A; • Q' where 



g-(i, ^), Q 

Notice that Q^,Q'^ are null vectors, Q^Qp 



^ (1- -p (18) 

Q'^Q'^ = 0. This is a remnant of the 
(5- function in Eq.(8) enforcing the mass-shell condition for the quarks. 

The ie's in the denominators in Eq.(15) will be let go to zero at this stage. The 
ie's can contribute to the imaginary part of the two-point function, corresponding 
to the Landau damping of the gluon field. The retarded rather than time-ordered 
propagators are appropriate for a discussion of damping effects; we shall obtain the 
correct imaginary part by an appropriate continuation later. For the moment, let 
us ignore the ie's. Using a, 13 from Eq.(5), we find, for the temperature-dependent 
part of r(^\ 



p(2) ^ _1 
^ 2 



dfi{k) 



d^q 1 1 
(27r)3 2pO 2qO 
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T{p,q) 



fco 



{rig 



Tip'.q') 



^P'pO 



qO + kO 



, T{p,q') 



, T{p',q) 
1"P + «9J + .0 + P 



(19) 



We see that the result is linear in the distribution functions, a property which 
holds in general for the n-point functions. Notice also how soft denominators like 
p° - (j° ± fc" arise from the time-integration. For \k\ small compared to the loop 
momentum |g|, we can write 



/-gO 



-k-Q p" -q'^ + k'^ ~k-Q' 

<o + 9° ± fc° ~ 2q° 



(20) 



T{p,q)~8q°\r{Ai-QA2-Q) 

T(p',g')-8/tr(Ai.Q'A2-g') (21) 

T{p\ q) ~ T{p, q') ~ V'tr(Ai • Q' A2 ■ Q + Ai ■ QA2 ■ Q' - 2Ai ■ A2) 

where Ai = A{k), A2 == A{k') and the remaining trace in the expressions for T's, 
denoted by tr, is over color indices. The difference of distribution functions can also 
be approximated as np- Ug c^ ^Q ■ k. (This is the extra factor of k/T mentioned as 
rule 5 above.) Using these results, Eq.(19) simplifies to 



r(2) = -i / d^i{k) 



_fq_ 

(2^ 



tr 



, dnAi-QA2-Q dn Ai ■ Q'A2 ■ Q\ ^^ - 
^dqO k-Q dqO k-Q' '^ 

—T^{Ai ■ Q'A2 ■Q + Ai-QA2-Q' - 2Ai • A2) 



(22) 



We have the result 



I d'<ipiQ) - - 1 d^ipiQ) 



(23) 



for any function / of Q, or Q'. We can further use 2Q -k = k- Q' - k- Q. Eq.(22) then 
simplifies to 



r(2) = -i / dfi{k) 



d^q 

(27r)s 



tr 



n + n 

r,0 



(2Ai ■QA2-Q-Ai- QA2 -Q'-Ai- Q' A2 • Q + 2Ai • A2) 



,J.^,.<3^,.0|_|_iL2^,.0M,-0'i| 



(24) 



The angular integration in Eq.(24) over the directions of q (or Q) can help in sim- 
plifying it further by virtue of 

I dfl (2Ai ■QA2-Q-Ai- QA2 -Q'-Ai- Q' A2 ■ Q + 2A1A2) ^ f dfl {2Ai ■ QA2 ■ Q') (25) 



Defining 



we can write Eq.(24) as 



A4 



AQ 

2 ' 



A_ 



A-Q' 



(26) 



r^^^--W..(fc)/|^^16tr 



k-Q' k ■ O 

Ai+A2-{n + n) - n^^Ai+A2+ - n^--^Ai_^2- 



(27) 



This expression becomes more transparent when written in coordinate space and 
finaUy with a Wick rotation to Euchdean space. Let us define the Green's functions 

G{xi,X2) ' 



P ■ Q (27r)4 

In terms of the nuU vectors Q = (1, q/q^), Q' = (1, -q/q^), we can introduce the hghtcone 
coordinates {u,v,x'^) as 

u = ^—, w=— ^, Q-x =0 (29) 

where Q = q/q°. We then have Q-d = du, Q'-d ^ dv We shall also introduce a Euclidean 
version of our results by the correspondence 

2u ^^ z, 2v ^^ z 

du^Q-d^ 2d^, dy^Q' -d^ 2dg (30) 

The Green's functions in Eq.(28) are the continuations of the Euclidean functions 

1 S(^^ix'[-x^) 



Ge{xi,X2) 



2-111 (zi - Z2) 



1 6i^){xJ-xl) 
G^ix.^x,)^— ^^^_^^^ (31) 

More precisely, the Green's functions in Eq.(28) obey the equations Q ■ dG{xi,x2) = 
-iS{xi ~ X2), Q' ■ dG'{xi,x2) = -i5{xi - X2). The Green's functions in Eq.(31) are 
the solutions to the corresponding Euclidean equations 2dzGE{xi - X2) = ~i5{xi - 
X2), 2dgG'^{xi - X2) = -i6{xi - X2)- This leads to the correspondence 

k-Q' 1 1 



k-Q TT Z12Z2I 

k-Q 1 1 

k-Q' TT Z12Z2I 

1 1 1 



k - Q 2111 zi2 

where z^ = Zi- Zj, etc. Eq.(27) for r^^^ can then be written as 






(32) 



(n + n) / d'^x A+{x)A^{x) 



rp(fzi(fz2A+{xi)A+{x2) _ f ,2 rj.(Pzi(fz2A^{xi)A^{x2) 

C ; n TT a X 

TT TT Z12Z21 J TT TT Z12Z21 



(33) 



This is finally starting to look nice. Let us define the functional I{A+) by the formula 

TT TT Zi2Z23---Znl 

'o tr( A , (T■.^A , (t^)) 

(34) 



I(A ) = iV-*-""^''" ^ -"^-^ ^'^^^ '^^^" t'^(^+(^i)-"^+(^«)) 



n J n TT Zi2Z23---Znl 

i f ,2 T <Pzi(Pz2 ti{A+{xi)A+{x2)) 

- dx ^^ + ••• 

2 J TT TX Z12Z2I 



We shall show later that I{A+) is related to the eikonal for a Chern-Simons theory. 
Eq.(33) for r^^^ can finally be written as 

d^q 1 



where 

K[A+,A-] = -16 



r(2) 



(n + n) I 4 



(27r)3 2q 



- X(2)[A+,yl_] 



d'^x tr (A+(x)A_(x)) +n iTTl{A+)+n i7r/(A_) 



(35a) 



(356) 



k'-'^''[A+,A^] in Eq.(35a) denotes terms in K which are quadratic in A. I is obtained 
from I hy z ^ z. 

Although we have introduced different distribution functions n, n for quarks and 
antiquarks, it is only n + n which is relevant at high temperatures. We see that, by 
virtue of /dO I{A+) = J dV, /(A_), we can write Eq.(35b) as 



K[A+,A^] = -16 



n + n 



d'^x tr{A+{x)A^{x)) +iTTl{A+)+iT:i{A^) 



(36) 



The integral over the magnitude of (fin Eqs.(35) can be easily carried out. With 
zero chemical potential, 

rp2 



r(2) 



127r 



dn 



d'^x tr(A+ AJ) +i7r/(2)(A+) +i7r/(2)(A_) 



(37) 



4. Properties of r[^] 

We can extract certain properties of T[A] from the power counting rules discussed 
earlier. There are two key properties which are important to our analysis. 

1) T[A] is gauge-invariant with respect to gauge transformations of the gauge 
potential A^ and is independent of the gauge-fixing used to define the gluon propa- 
gators ^'^°. 

We can understand how the kinematics of hard thermal loops can lead to gauge 
invariance. The thermal propagators satisfy the same differential equations as the 
zero temperature propagators. This is evident from Eqs.(8,9b); only the choice 
of the homogeneous solution, equivalently boundary condition, is different. As a 
result, the generating functional of one-particle irreducible vertices, viz., T[A,c,c,q,q\ 
obeys the standard BRST Ward identities, c and c are the ghost and antighost 
fields respectively. For the usual gauge fixing term \{d ■ AY, we thus have ^^ 



d\ 



5T 6V 

Ja^Ik^ 



srsr ^^^ ^^sr 



(38) 



Kfj^ and L are sources, corresponding to the terms K'^{Df^c)°-, f^^'^L'^c^c" added to 
the action. In the hard thermal loop approximation, terms involving the ghosts 
are subdominant. Recall that the ghost-gluon coupling involves /"'"=A°^(9^c'')c'=; the 
derivative is on the antighost field. In a diagram with external ghosts, this is a power 
of external momentum and therefore such a diagram is smaller compared to a similar 
diagram with the ghost lines replaced by gluons where there are contributions with 
all derivatives on the internal lines which give powers of T. (See figs. (3).) The 
terms ^, |^ in Eq.(38) are thus negligible in the hard thermal loop approximation. 
Further, we are interested only in one-loop terms. The identity (38) then gives 
the gauge invariance of r. Thus effectively, in a high-T-expansion, the leading 
term T[A] which is proportional to T^ is gauge-invariant. Notice also that since the 
thermal contribution to the propagator is on-shell, the T-dependent part of a one- 
loop diagram is classical and so it is not surprising that the BRST Ward identities 
reduce to the statement of gauge invariance. The fact that r does not depend on the 
gauge choice for the gluon propagators can be seen by similar arguments; identities 
for the variation of r for changes in A can be written down and simplified ^'^°. 
2) T[A] has the form ^^ 



■2 



r = (iv + iAr.)^ 



d'^x 27rAgAS + I d^ W{A ■ Q) 



(39) 



Here Q^ is the null vector (1,Q) with Q^ = 1, Q^Q^ = 0. The AgAg-term is the lowest 
order Debye screening effect; it is clearly an electrostatic mass term. The key point 
about Eq.(39) is that, for each Q, T[A] involves essentially only two components of 
the gauge potential, Aq and A ■ Q. The dr^-integration will bring in all components 
of the potential, but for each Q only two components are needed. We have already 
seen this structure in the explicit calculation of the two-point function, where Q was 
the angular part of the loop-momentum q. The dfi-integration in Eq.(39) is over 
the orientations of Q and is the unfinished part of the loop integration, after the 
integration over the modulus of the loop momentum has been done. The structure 
of Eq.(39) can be seen by analysis of diagrams again. For example, for the diagram 
of fig. (3b), with the derivative gluon coupling, since p ■ A ^ q ■ A, the possible tensor 
structures are q'^A'^ and {q-A)^. The former is zero since the propagator involves the 
(5-function 6{q^). Writing q = \q\{l,Q) and carrying out the |q|-integration, we are left 
with a structure like Eq.(39). This argument generalizes to diagrams with arbitrary 
number of external gluons. 

Given these two properties of r one can determine W and hence r simply by the 
requirement of gauge invariance ^^. The condition for gauge invariance of r is 

dn 6W = Att I d'^x A^uj" (40) 

where Ag is the time-derivative of A^ and oj = -it°-uj°- is the parameter of the gauge 
transformation, i.e., 5A^ = d^uj + [Ap,cj]. Eq.(40) can be realized by 

5W ^ f d^x A" -Qu"-. (41) 



One can check that Eq.(41) is indeed the way gauge invariance is reahzed by anal- 
ysis of the diagrams. It is clearly so for the two-point function from our explicit 
calculation. We now rewrite Eq.(41), using 

f , SW SW 

SW^^J dHQ • d^^^^ + [A . Q, j^^^]ru.^ (42) 

as 

where 

SW 

We have used the lightcone coordinates from Eq.(29). We now make the Wick 
rotation as in Eq.(30) which gives A+ = ^A- Q -^ A^; also rename / as 

The condition of gauge invariance, Eq.(43), then becomes 

d,A, - d,a, + [a,, A,] = 0. (46) 

If ^2, flj are thought of as the gauge potentials of another gauge theory, we 
see that Eq.(46) is the vanishing of the field strength or curvature F^z- The gauge 
theory whose equations of motion say that the field strengths vanish is the Chern- 
Simons theory. We shall therefore put aside Eq.(46) for a moment and turn to a 
short digression on the Chern-Simons theory, returning to Eq.(46) and its solution 
in section 6. Of course, an understanding of Chern-Simons theory is not absolutely 
essential to solving Eq.(46). One can simply solve Eq.(46) and regard Chern-Simons 
theory as an interpretation of the mathematical steps along the way. However 
Chern-Simons theory does illuminate many of the nice geometrical properties of 
the final result and is a worthwhile digression. 

5. Chern-Simons and WZNW Theories 

The Chern-Simons theory is a gauge theory in two space (and one time) dimen- 
sions ^^'^'*. The action is given by 



lMx[ti,tj] 



5* - ^ / d^x e^"^" tr(a^9^a„ + |a^a^a„). (47) 



Here a^ is the Lie algebra valued gauge potential, a^ = -if^a"^. k is a constant whose 
precise value we do not need to specify at this stage. We shall consider the spatial 
manifold to be R^, or C since we shall be using complex coordinates z ^ x + iy, 
z = X- iy. (Actually, we have sufficient regularity conditions at spatial infinity that 



we may take M to be the Riemann sphere.) The equations of motion for the theory 
are 

Ffiy = dfj_a^ - dya^j_ + [a^, a^] = 0. (48) 

The theory is best analyzed, for our purposes, in the gauge where a^ is set to 
zero. In this gauge, the equations of motion (48) teU us that a^, oj are independent 
of time, but must satisfy the constraint 

Fj^ = dsaz - dzttz + [a-,, az] = 0. (49) 

This constraint is just the Gauss law of the CS gauge theory. It can be solved for 
flj as a function of a^, at least as a power series in a^. The result is 

J n TT (Z - Zi)(Zi - Z2) ■ ■ .{Zn- Z) 

This can be easily checked using dz{j^) = nS'-^^z - z'). 
In the flo = gauge, the action becomes 



dtd X ir^azdoaz). (51) 

This shows that a^ is essentially canonically conjugate to az. In fact in carrying out 
a variation of S*, we find the surface term 6{tf) - 9{ti), where 

e=— I (fx ir{az5az). (52) 

'^ Jm 

(We assume agSuz to vanish at spatial infinity.) 9 is the canonical one- form of the 
CS theory. (This is so by definition; the canonical one-form in any theory can 
be defined by 9, where the surface term in the variation of the action is 9f - 9i, 
the subscripts referring to the final and initial data surfaces.) 9 is the analogue 
of ftdx* of point-particle mechanics; for an action S == J dtdx [^^ - V{x)], we would 
find 9 = mxi6x^ ^ piSx\ We can make another variation of 9, antisymmetrized with 
respect to the variation Saz, denoted by the wedge product sign, and write 

Lo ^. 59 = — / d X tr{5az A 5az) 
TT Jm 



d^xd^x' ljab{x, x') S^'^ix) A 6^^{x') (53) 

where 

, „ in f S{x-x')S''''\ ,^,. 

u:Mx,x)^--^_^^^_^,^^,, ) ^^^) 

and 6£,^ = {5a%,5al). uj defined by Eqs. (53,54) is called the symplectic structure and 
is of course the analogue of dpidx^ of particle mechanics. The inverse of ujab gives 
the Poisson brackets, the commutators being i times the Poisson brackets. For our 
case we get 

[e^(x),^^(x')]=z(^-^)-^^(x,a;') 



or 

[a%{x),a\{x)]^-5'^'5^^Hx-x'). (55) 

K 

One does not have to go through 9 and w to arrive at Eq.(55). One could simply use 
the fact that, from Eq.(51) the canonical momenta are tt = oj, tt = 0. This is thus a 
constrained system in the Dirac sense and using the theory of constraints one can 
derive Eq.(55). The procedure of using ujab{x,x') is quicker. 

In the expression (52) for 9, ag is independent of a^. We can however express aj as 
a function of az via the constraint Eq.(49) or equivalently Eq.(50) and functionally 
integrate 9. In other words, we define I{az) by 



51 — — / d X tT:[az{az)5az]- (56) 

71" J 

The solution for / is given by 

j_ . V^ (-1)" [ <Pzi (fzn iT{az{zi,Zi) . . .az{Zn,Zn)) , . 

^ n J TT TT Z12Z23 ■ ■ ■ Zn-lnZnl 

The quantity / has a rather simple interpretation. For one-dimensional point- 
particle mechanics, 9, as we mentioned earlier, is given by pdx. p is independent of 
X to begin with, but we can express it as a function of x via a constraint such as of 
fixed energy, e.g., ^ + V{x) = E. Integral of 6* = pdx then gives Hamilton's principal 
function or the eikonal, familiar as the exponent for the WKB wave functions of 
one-dimensional quantum mechanics. We have an analogous situation with Eq.(49) 
expressing ag as a function of a^. / is thus an eikonal of the CS theory ^'''^^. 

/ is in fact the Wess-Zumino-Novikov-Witten (WZNW) action ^^. We can write 
the gauge potential Uz as a^ = -dzUU~^ where U is in general not unitary; it is an 
SL{N, C) matrix for gauge group SU{N). Notice that since dz has an inverse by virtue 
of dz (z-z') = T^^^^^{z - z'), such a U can be constructed for any az, at least as a power 
series in a^. / can then be written as / = -iK.SwzNw^) where 

SwzNw{V) = ^ f d^x tr{dzUdzU-') - -^ / d'^x e^'"'tr{U-'d^UU-'d,UU-'do.U). (58) 

The second term, the so-called Wess-Zumino (WZ) term, involves an extension 
of U into a three-dimensional space. We take M^ = m x [0,1] with U{z,z,0) = 1, 
U{z,z,l) = U{z,z). The integrand of the WZ term becomes a total derivative once a 
(local) parametrization is chosen for U. This ensures that the physics is independent 
of small changes in how one makes the extension of U into the extra dimension. If we 
have two globally inequivalent extensions, say Ui and C/2, the difference SwzNwiUi)- 
SwzNw{U2) = '2TrQ[g] where Q[g] is the winding number of the map g : S^ —^ G; G is the 
space in which the matrices U take values, SL{N, C) in the present case and g ^ Ui 
on the upper hemisphere of S^, which can be taken as one copy of M^ and g = C/2 on 
the lower hemisphere, taken as another copy of M^. Since Q is an integer, e' 
is independent of globally different extensions if k is an integer ^^. 



The WZNW-action also obeys the foUowing property, sometimes caUed the 
Polyakov-Wiegmann formula ^^. 



S{U) - - I Tr{h-^d,h d,U U-^) 



S{hU) = S{h) + S{U) / TT:{h-^dzh d,U U-^) (59) 



The crucial point here is that in the cross-term h has only z-derivative and U has 
only z-derivative. (They become 9_, d+ in Minkowski space.) 

The relationship of Swznw to / is easily seen by considering variations. Under 
the variation U -^ e'^C/ ~ (1 + tp)U, we find 5a^_ = -D^^ip = -{dzip + [a^, tp]) and 



SS- 



WZNW 



- I (fx tr{diif>a^). (60) 



Partially integrating and using F^z = and Sa^ = -D^ip, we find that Swznw obeys 
Eq.(56) except for a factor (-^k), thus identifying / = -ikSwznw- 
Another quantity of interest is 

1 



K 



K d X ti^azaz) + inl^az) + iTTl{ag) 



(61) 



It is easily checked that this is gauge-invariant. K has a nice interpretation as a 
Kahler potential, which I shall not go into here ^. 

Finally notice that the eikonal / of Eq.(57) may be regarded as the expansion 
in powers of a^ of the logarithm of the functional determinant of Dz = dz + az] 
i.e., / = {-iK)\ogdetDz = -in TilogDz- The expression K of Eq.(61) is then given by 
-KTT\og{DzDz). The expansion of this expression in powers of the potential obviously 
gives / and /. The extra term / ^tr{azaz) is precisely the local counterterm needed 
to give a gauge invariantly regulated meaning to TT\og{DzDz). 

6. Solution for r[^] 

We now return to Eqs. (45,46) for the quark-gluon plasma. We can rewrite 
Eq.(45) as 

5W ^A I ct^x ty:{az5Az) - 5 j ct^x A^A^. (62) 

Comparing with Eq.(56), we see that, since a^, Az obey the constraint Eq.(46), the 
solution is related to the eikonal /. The difference here is that Az and a^ depend on 
all four coordinates x^, not just z, z. However, there are no derivatives with respect 
to the transverse coordinates x'^ in Eq.(46) and hence the solution for a^ in terms 
of Qz is the same as in Eq.(50), with a^ depending on x^ in addition to z, z. The 
x^-argument of all az factors is the same, i.e., 

E. ,^„_l [ (Pzi CpZn Az{zi, Zi, x'^)Az{z2, Z2, X^) . . . AzjZn, Z„, x'^) 

J n TT (Z - Zi)Zi2 . . .Zn-ln[Zn- Z) 

For the eikonal we get the same expression as Eq.(57) but with integration over the 
transverse coordinates, i.e., 

j_.^.^y (-1)" f^2^T d^z^ d^z^ tT{Azix^)...Az{xr,)) ^ ^g^^ 

^ n J TT " ' TT Z12Z23 . . . Z„l 



Since it is not relevant to the present discussion, we have, for the moment, set k = 1. 
We then find the solution to Eq.(62) as 



W = -AmI{A^) - 



d'^x At At. 



(65) 



(Strictly speaking a^ and A^ are not complex conjugates; the analogy holds better 
with a Chern-Simons theory of complex gauge group. However, we are only using 
the Chern-Simons analogy to obtain Eq.(65). The expression for r[^] can also be 
directly checked to be a solution to Eqs.(45,46).) 
From Eq.(39), r is given by 



r = (iv + ^Nf) 



127r 



d-^x 27rAgAg - dn{ d^x A^Al + AmI{A;)} 



We can now use the identity 



I dn tv{A,A-,) = -\ 



2t,AIAI 



dn {AtAD 



which follows by straightforward dfi-integration, to write 



-{N+\Nf) 



y2 



2- ^'6^ 
2 



dn 



d'^x tr{A^Ag) + inl^A^) + ini^Ag) 



{N+\Nf)^ j dnK{A,,A,) 



(66) 

(67) 

(68a) 
(686) 



where K is given by Eq.(61) with k = 1 and with the additional integration over 
coordinates x^ transverse to Q. If we write A^ = -dzUU~^ and Ag = W~ dgU, we can 
write Eq.(68) in terms of Swznw, using Eq.(59), as 



-{N+^Nf)^Swznw{U^U). 



(69) 



(Again, a suitable additional integration over the transverse coordinates is under- 
stood.) 

From Eqs. (41,42), it may seem that our solution is ambiguous up to the addi- 
tion of a purely gauge invariant term. But for hard thermal loops, the additional 
structure that W depends only on A^ tells us that the gauge invariant piece must 
obey Dzjj^ = 0. Since d^ is invertible, at least perturbatively, there is no nontrivial 
solution to this equation. Thus Eq.(66) or Eq.(68) is the unique solution and r so 
defined must indeed be the generator of hard thermal loops. 

In the last section, we noted that K{Az,Ag) can be considered as TAogiD^Dg). 
Since D^ and D-, are the chiral Dirac operators in two dimensions, TTlog{DzDz) is 
clearly the photon mass term of the Schwinger model, for Abelian gauge fields. 
More generally, for non-Abelian fields as well, we can consider Trlog(Z)zZ)j) as a 
gauge-invariant mass term. It is perhaps fitting that the gauge-invariant Debye 
screening mass term in four-dimensional QCD is given by suitable integrations of 



such a two-dimensional mass term (with, of course, the additional x^-dependence) . 
Equally appropriately, the Chern-Simons term, from which we derive K{Az,Ag), is 
also the mass term for gauge fields in three dimensions. 

We close this section with two remarks on r. The CS theory of Eq.(47) violates 
parity. We see that this parity violation disappears, as indeed it should, by integra- 
tion over the orientations of Q, for the QCD case. Alternatively, expression (68a) 
is manifestly parity-symmetric with A^ ^ Ag (or A+ ^ A_ in Minkowski space), 
Q ^ Q' under parity. Secondly, instead of setting k = 1 and then having a prefac- 
tor (N + \Nf)^ in Eq.(68), we could simply choose k = (iV + \Nf)^. For WZNW 
actions, if U has a non-Abelian unitary part, any action we use must be an inte- 
ger times SwzNw, as we have already seen. In our case, r involves SwzNwiV^U) or 
SwzNw{H) where H is hermitian. The third homotopy group is trivial for the space 
of hermitian matrices and so, as expected, there is no argument for quantization of 
the coefficient, which is {N + \Nf)^ for us. 

7. Continuation to Minkowski space 

We shall now consider how the above results can be continued back to Minkowski 
space ^. Evidently, we expect A^ ^ A+, A^ —> A^ and d^ —> d+, dg —> d-. The inverses 
of 5±, however, require specification of boundary conditions or ie-prescriptions. This 
can be easily seen in momentum space where the inverses have denominators like 
k ■ Q which can vanish; one must give a prescription on how these singularities are to 
be handled for the ^"-integration. The appropriate ie-prescription depends on the 
physical context and the quantity being considered. Very often, one is interested 
in the evolution of fields in a plasma, which can be described by the operator field 
equations. For this case, the retarded Green's functions are the appropriate ones. 
Let me illustrate this with a simple example, viz., a plasma in electrodynamics. 
The field equations can be written as ^® 

92 A^ = is-^^ = jt' (70) 

oA'iJ' 

(We have chosen the Feynman gauge for simplicity.) S is the scattering operator 
considered as a function of the incoming field A™ . We can expand the current as a 
function of the field A^\ obtaining to linear order in the field, 

J" = .r -^jd^y e(a;° - y^'WixlfiyW^^iy) (71) 

where j^ is the current in the interaction picture, say qi^^qi in terms of the fernion 
field qi in the interaction picture. Notice that the we have the retarded commutator 
of the currents j^. Eq.(71) is an operator equation which is generally true. We 
can take any kind of expectation value of this equation; the thermal result which is 
the induced current in the plasma is obtained by taking a thermal average. (The 
result is the Kubo formula.) To the order we are interested in, the incoming and 



interacting fields can be considered the same and Eq.(70) gives its evolution. The 
appropriate boundary condition for the induced current which governs the evolution 
of the fields, we see, is the retarded condition. Of course, the induced current can 
also be calculated in terms of the induced or effective action r as -j^-- The strategy 
for continuation to Minkowski space is thus to calculate the induced current in 
Euclidean space and then to use ies appropriately in denominators involving fc • Q's 
so as to get the retarded condition. If J^ is expanded to higher orders in the field, we 
encounter multiple retarded commutators. The ie's must be inserted appropriately 
so that we get this retardation structure. With this understanding, we can write 
the equations for the evolution of fields of soft momenta in the quark-gluon plasma 
as 

D^F"'''' = J^'" (72a) 



oo „ 
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jrik) = - dn 
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here 



Tr ( (^ii^)A_(fcO + A+{k,)iA^) 1 J„.i 



+ i -(2*)"-iTr ( (-!^)A+(fci) • • • A+iK) ] F{h, ■ ■ ■ fc„) + (Q ^ Q') 



(72c) 



Fih,---kn)=J2 



Qi 
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{qt - qa){qi ~ qi)--- (qt - qt-i){qi - qi+i) • • • (<?» - qn) 









We may reexpress the current as 

J"^" = -^ /" dl^ Tr [(-te) {{a+ - A+)Q"' + {Q' ^ Q)}] 

where a+ is defined by 

d-a+ - d+A^ + [A_,a+] = 



(73a) 
(735) 



(74) 



(75) 



The retarded condition is to be used in solving Eq.(75). 

It is easy to check that the imaginary part of the current given by these equations 
gives the Landau damping effects. To linear order in the field AJ^, we have J^ = 

Jd^yU-iix,y)A--iy),with 



u;lix,y)^6^''J^ e^^^^-y^ n^.(fc) 



n^,.(fc) = - 



27r 



47rg^offi^o - fco I dO. 



k ■ Q ^ ie 



(76) 



The imaginary part exists only for spacelike momenta and is evidently given by 
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(77) 



' - |fc|2 

As I said before, to the accuracy we are interested in here, the distinction be- 
tween the interacting and incoming field is irrelevant. One can, of course, go beyond 
this and include higher order corrections, using the time-contour approach due to 
Schwinger and Bakshi and Mahanthappa, which was rediscovered by Keldysh ^^. 
The operator approach is conceptually clearer in indicating why the retarded con- 
dition is appropriate, but a functional rewriting is useful in going to higher orders. 

We define 

Trp Tc cxp{iSint + iA-ri) 

^[^1 = TV ^^^^ 

where the time-integral goes from -oo to oo, folds back and goes from oo to -oo; Tc 
denotes ordering along this time-contour, p is the thermal density matrix and ry^ is 
a source function. One can represent Z[r]] as a functional integral 

Z[t]] = / dfiiA, c, c) exp{iSc{A, c, c) + i^^r;^) (79) 

The action is again defined on the time-contour. The Green's functions which arise 
in perturbatively integrating out the fields are to be taken as the real-time thermal 
Green's functions, as in Eqs.(4,9). In this time-contour version, the field is given by 
is 

—i oZ 

<^M>=^T TT 80 

^ Z Sr]t^{x) ^ ' 

We differentiate with respect to r\ on the first branch for the time-contour. It is 
evident that we can now define a generator for the one-particle-irreducible graphs 
rc[A], such that ^^ = -77. The evolution of a field configuration in the medium 
is given by ^ == 0; we can now set the source to zero. Separating out the term 
\{di,Ay)'^, which gives 9^A^, we get the current as 

< J" >= ^ (81) 

Vc = J(j ^(d^Afj^)^ + Tp. Eq.(81) is the same as the expectation value of our operator 
definition, with one-particle reducible terms summed up and expressed in terms of 
< A"^ > rather than AJ^™. In writing out Eq.(81), we encounter A's on the second 
branch of time. From the operator definition, we see that fields on the second branch 



are the conjugates of fields on the first branch. (In the usual time-contour approach, 
one considers equations for the two-point functions, here we consider equations for 
the fields; this is the only difference.) 

To summarize, the functional transcription of the operator analysis gives the 
following. The evolution of fields in the medium is given by ^ = 0, with fields on 
the second branch being conjugates of fields on the first branch. One can check that 
this gives the retarded prescription we use. For example, upto linear order in the 
field Ai^, we get. 



Jc 
= -M / < T{j^{x)My)) >+ < My)j^{x) > A'^iy) , (82) 

= -tjeix'^ - yO) < lUx),jAy)] > A-'iy) 
which agrees with Eq.(71). 

8. Derivation of r[A] from the Boltzmann equation 

The equations of motion we obtained and the expression for the induced current 
have also been derived using suitable truncations of the Schwinger-Dyson equations 
^", from effective actions for composite operators ^^ and from kinetic theory ^'^. 
These are all related approaches; here I shall briefiy indicate the kinetic theory 
approach. The classical equations of motion for non-Abelian particles, the so-called 
Wong equations ^^, are 

m— — = pf" (83a 

(IT 

^-^ = m^F^P. (836) 

"^^ - -.9/"''(p^4)9' (83c) 

Here q°- represents the classical color charge of the particle. (The color degrees of 
freedom can also be described in a phase space way; the appropriate space is the 
Lie group modulo the maximal torus '^^.) 

For particles obeying the Wong equations, the coUisionless Boltzmann equation 
for the distribution function f{x,p,q) is given by ^^ 



pM 



^^-''^"^^^^wr'-^^'^^'^i: 



f{x,p,q)=0 (84) 



The Boltzmann equation is invariant under gauge transformations g^^ —> UgAf^U^^- 
Ud^U^^, f{x,p,q) —> f{x,p,UqU^^). We then seek a perturbative solution of the form 
/ = /'-"^ + ff/^^' + •••, where /("' = Up is the equilibrium distribution, appropriately 
chosen for bosons and fermions. The Boltzmann equation, to the first order, gives 



P^ 






f^'^^p'^qaF^^^J^"^ (85) 



The color current, to leading order, is 

J^^g^ I [dq] Pf.q'^f'^ (86) 



[dq] is a measure for integration over the color charges which can be constructed in 
terms of the group coordinates modulo the maximal torus. By virtue of Eq.(85), J° 
satisfies the equation 



(p ■ DJ'^r - sVp'^^L^ (I dq g«gb/(")) (87) 

We now integrate this equation over the magnitude of p. Defining 



j^ = dn J 
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we see that Eq.(87) for the current becomes the zero- curvature condition Eq.(75); 
the rest of the analysis is as before, leading to Eqs. (72,73). This derivation from 
kinetic theory shows that the hard thermal loops are classical; they include thermal 
fiuctuations, not quantum fiuctuations. Of course, we have already seen this from 
another point of view, viz., the thermal part of the propagator is on-shell and so, 
the T-dependent part of a one-loop diagram describes the tree-level absorption and 
emission of particles from the heat bath. 

9. An auxiliary field and Hamiltonian analysis 

We now turn to the question of how we can use the effective action T[A]. There 
are two related but different contexts in which we need the expression for r[A]. The 
first is in setting up thermal perturbation theory. The version of T[A] as given in 
Eqs. (64,68) is probably best suited for this purpose. The resummed perturbation 
theory can be set up as follows. We introduce a splitting of the Yang-Mills action 
as 

S = So — cT 
So^ Jd^x ( -iF2) +r (89) 

We define propagators and vertices and start off the perturbative expansion using 
So- cT will be treated as a 'counterterm', nominally one order higher in the thermal 
loops than So- Eventually of course c is taken to be 1, so that we are only achieving 
a rearrangement of terms in the perturbative expansion. (As usual we must have 
gauge fixing and ghost terms.) Using this procedure one can calculate quantities 



which require resummations such as the gluon decay rate in the plasma ^. In 
Eq.(89) we have not displayed the quark terms. In addition to the quark kinetic 
energy terms, there are hard thermal loops which give a T-dependent mass to the 
quarks ^^^^^. We have not discussed this so far, since Chern-Simons theory does not 
give any new insights on this question. This mass term is actually given by 



r(«)[g,A] = / d'^xd'^y q{x)Y.{x,y)q{y) 

n.,y)^'-^ldn,.QF{x,y) (90) 

{iQ-D)F{x,y) = S^^\x-y) 

where Cp = f^f^ is the value of the quadratic Casimir for the quark representation. 
This action must also be added to and subtracted from the QCD action and treated 
in a manner analogous to r[A] to set up the perturbation theory properly. 

In addition to its role in rearrangement of perturbation theory, we can also use 
r[A] added to the usual Yang-Mills action, viz., 

Seff = J d^x i-lF^^Fr) + r[A], (91) 

as an effective action for the soft modes. (This is the spirit of Eqs.(72,73).) The 
nonlocality of r[A] makes it somewhat difficult to handle in this context. It is useful 
to rewrite T[A] using an auxiliary field which makes the equations of motion local ^. 
The auxiliary field is also useful in setting up a Hamitonian analysis of the dynamics 
of the soft modes. The auxiliary field we use will be an SU{N)-uiatnx field G{x, Q) 
which is a function of x and Q, i.e., defined on M'^ x 5*^, M'^ being Minkowski space. 
Further G{x,Q) must satisfy the condition G^x,Q) = G{x,-Q). The action is given 
by 

S = f -^F^ + K f dfl d^x^ SwzNwiG) + - I d'^x T:i{G-^d-G A+ 



A^d+G G-^ + A+G-^A^G - A+A^ 



(92) 



where Swznw{G) is the WZNW action of Eq.(58) with U replaced by G. The quantity 
in the square brackets in Eq.(92) is the gauged WZNW action ^'^. It is invariant 
under gauge transformations with G transforming as G ^ G" = h{x)G h^^ix), h{x) e 
SU{N). The equations of motion for the action (92) are 

d+A--d-a+ + [a+,A-]=0 (93a) 

a+ = GA+G-^ - d+G G"^ (936) 

{Df.Fi"')'' -r" =0 (94a) 



r'' = -^ fdn TT[{{-ima+ - A+)Q"'} + (Q' ^ Q)] (946) 

Ztt J 

= -^ fdn Tr[(-zt"){G-iZ?_G Q" - D+G G-^Q'"}] (94c) 

27r J 

Clearly these equations are equivalent to Eqs. (72,74,75); the only difference is that 
the equation defining a+ in Eq.(74), viz. Eq.(75), is now obtained as the equation of 
motion (93a) for G. Notice that the current in Eq.(94c) looks like the current of a 
matter field; thus except for the fact that G depends on Q, QCD, with hard thermal 
loops added, is no stranger than Yang-Mills theory coupled to a matter field. The 
equations of motion for G has no independent solutions, so that we are not changing 
the physical degrees of freedom by introducing the auxiliary field. This can be seen 
as follows. We can parametrize the potentials A± in terms an S't/(Af)-matrix V(x, Q) 
as 

A+ = -d+V y-^, A- = -d-V V'-^ (95) 

where V'{x,Q) = V{x,-Q). The general solution to Eq.(93) is then given by 

G{x,Q) = V{x,-Q)B{x+,x'^,Q)C{x~,x'^,Q)V-\x,Q) (96) 

where G is an arbitrary 5'C/(iV)-matrix depending on the variables indicated and B 
is given by C with Q —> -Q. The matrices B, C represent the new or independent 
degrees of freedom for the field G. Notice, however, that the parametrization (95) of 
the potentials has redundant variables; a transformation V —> VU{x^,x'^,Q), with a 
corresponding change in V', leaves the potentials invariant. Further, ordinary gauge 
transformations act on the matrix V" as y ^ h{x)V. Thus the physical subspace for 
the components A± is given in terms of the matrices V with the identifications 

V{x,Q)^h{x)V{x,Q)U{x-,x'^,Q) (97) 

The gauge freedom of multiplying V"'s by matrices which do not depend on x+, viz., 
C/'s in Eq.(97), shows that the we can reduce G to just V{x,Q)V~^{x,Q). There are 
thus no new real dynamical degrees of freedom in G. The action can be simplified 
for this G to 

r[A,G] = -K [ dn d^x'^ SwzNw{V-\x,-Q)V{x,Q)) (98) 



This agrees with the expression obtained by Feynman graph evaluation of the hard 
thermal loops. 

The Hamiltonian for the effective action (92) can be obtained by following the 
usual analysis for the WZNW-action; it is given by ^ 

n= f d^^ \ ^^ + ^ f dnTr\^{DoG DoG-^) + {Q-DG Q-DG-^^-A^gA (99) 

where F^^ = E^, F^^ = e.^kBl and 

g" = {D-EY + -^ I dmr[{-ie){G-^D_G - D+G G-^)] (100) 



tj'" = is the Gauss law of the theory; it is also the time component of the equation 
of motion for the gauge field. Expression (99) makes it clear that the Hamiltonian is 
positive for all configurations which are physical, i.e., obey the Gauss law. (It may 
be worth recalling that, even for the Maxwell theory, the canonical Hamiltonian is 
positive only for physical configurations.) We find, comfortingly, that the effective 
theory for the soft modes has positive energy. 

The Hamiltonian analysis is most easily carried out, as in the usual cases, in the 
gauge Aq = 0. We start by defining the currents 

J+ = ^D+G G-^ = {-ie)Jl 

47r \ ' + 

J- = -^G-^D^G = {-ie^J" (101) 

By virtue of the property G^^{x,Q) = G{x,-Q), these are related by J+{x,-Q) = 
J-{x,Q). The Hamiltonian can be written in terms of these currents as 

H= I dPx 1 ^ ^^ +— f dn (J^Jl + J1J1)\ (102) 

We have chosen the Af, ^ gauge; Gauss law must henceforth be imposed as a 
constraint. ( For fixed Q, the integrand of the second term involving the square of 
J+ and the square of J_ is the Sugawara form of the Hamiltonian, well-known in 
the context of two-dimensional current algebras.) The equations of motion in the 
Aq = gauge are 

Ef ^dnAf (103a) 

doEf + e,,k{DjBkT = j d^ Q,{Jl - r) (1036) 

Equation (103a) is just the definition of the electric field; however, in a Hamiltonian 
approach, it is an equation of motion and we have displayed it as such. 

The commutation rules must be such that Eqs.(103) follow as the Heisenberg 
equations of motion for the Hamiltonian (102). Knowing the current algebra of 
the WZNW-theory, we can make a guess as to what the appropriate commutation 
rules are for our problem and verify them by checking that they lead to Eqs.(103) 
starting from the Hamiltonian (102). The commutation rules are then seen to be 

[E^{x),A){x')] =-i5''^5ij5{x-x') (104a) 

[ El{x), Jlix') ] = ±i^Q^5''^5{x-x') (1046) 

47r 

[ jj(x,g),4(f' ,g') ] ^ir''^Ji5{x~x')5{QM') 

T^Qr{D,)t5{x-x')5{Q,Q') (104c) 

47r 

[ J°(x,Q),J^(x',Q') ] =0 (104d) 



All other commutators vanish. S{Q, Q') stands for the i5-function on the sphere 
corresponding to the unit vector Q, i.e., j dQ.Qi 5{Q,Q')f{Q') = f{Q). 

The algebra (104) is an interesting extension of the usual WZNW-current alge- 
bra; we are in four dimensions and have a gauge field as well. We can check that 
the commutation rules (104) obey the Jacobi identity. Since we are postulating 
the algebra, this is a necessary check. Of course, commutation rules can also be 
obtained from the action by standard quantization procedures ^. The condition 
J+{x - Q) - Jl{x,Q) = has to be imposed as a constraint, just like the Gauss law. 

10. Plasma waves 

Long wavelength and low frequency plasma waves are the classical solutions of 
the effective theory Eq.(91). It also includes effects such as the screening of Coulomb 
fields. These features can be seen by examining the Abelian case or electrodynamics. 
In this case, the terms in r which are cubic or higher order in A^ are zero and in 
terms of the Fourier components of A^, we can write 

Seff - i / 7^ A^{~k)M^^{k)A,{k) (105a) 

where 

i^ r f /OM/o^ 

(1056) 

k'^Mfj.i, = in accordance with the requirement of gauge invariance. k = Npe'^T'^/6 for 
electrodynamics; we have restored the coupling constant e in this expression. For 
the plasma waves, we will see that the wavevector is timelike; hence the imaginary 
part of Mfj^i, is irrelevant. 

We can split ^^ into a gauge dependent part and gauge invariant components 
as 

A^ = fc^A(fc) + a^ + l3^ (106) 

where A shifts under gauge transformations and a^, /3^ are gauge invariant. We take 

, k^ ,j, ko (3) . ^2 



1 [ d'k 


^i~k)M^''{k)A,ik) 


k'^kn + f 

ZTT 


[47rg^V" k'Jdn f^^j 



h-l Z^2 ^ / -To Z^2 Z^2 

/3o = A = ef^aA, A = 1,2. (107) 

(Here we are considering fields off-shell and so fc° ^ \k\.) The Ci's form a triad of 
spatial unit vectors which may be taken as 



s« = (e.,4f'0)' e<^^-(^i,-J^), ^ = 1,2, (108) 

V4 Jfc|fc2 




where fc|, = kf 



1^2 • 



Notice that kici ' = 0, A = 1,2. (/) and a\ are the gauge invariant 



degrees of freedom in Eq.(106). When the mode decomposition (106) is used in 
Eq.(105) we get 



^eff - 



where 



(2^)4 



A(-fc)(- '' 



fc2 



^)M^(fc)/3,(fc) 



(f){-k)M^{k)(f>{k) 



M'^{k) = kl-k^ -n^{ka,k) 



i L'Y> hi 



F B 2\k\ 



(pj^+p^r (109) 



M^{k) 



(ki-k^) 






rtl 






fc2 




L = log( 



ko 



(110) 



(111) 



'ko-\k\ 

We have also included an interaction term with a conserved source J^ in Eq.(109); 
i.e., we include j A^J^ and simplify it using Eq.(106). From Eq.(108) we see that the 
interaction between charges in the plasma is governed by {M^y^ which shows the 
Debye screening with a Debye mass m^ = \/2k. The action (109) can also give free 
wavelike solutions. The dispersion rules for these plasma waves would be M'^ = 
for the transverse waves and ( '' Z^" )M^ = for the longitudinal waves ^s. (The extra 



fc2 



factor multiplying M^ is from rewriting 4>{k) in terms of the potential.) 

Non-Abelian plasma waves can be defined in a similar way as propagating so- 
lutions to the equations of motion given by Eq.(91). For plasma waves in Abelian 
subgroups, only terms in r upto the quadratic order in A^ are important; the anal- 
ysis is the same as for the Abelian plasma waves with m|, = 2k = (A^ + ^Nf)^-^ and 
K = (TV + ^Nf)^ in M^{k). 

It is interesting to consider these waves also from the Hamiltonian point of view. 
We can choose Ao = gauge and introduce the parametrization 



A, 



(3) , (A) 

e,- 'a + e, 'ax 



E..^^Ii- 



ef^n. 



We have dropped the color indices, since we are in an Abelian subgroup. 
Hamiltonian becomes 



(112a) 
(1126) 

The 



n= (fx 



tf + n^ + djaxdjax + — I df^(J+ + J't) 



(113) 



The commutation rules are 



[n(x),a(27)] = -i5{x-y), 
[^(f),J±(y)]=±^|^g•e(3)5(f-y), 



[\lx{S).a'M\ - -^5xx'5{x-y) (114a) 

[Tlx{x),J±{y)]^±i^Q-e^^H{x-y) (1146) 



The creation operators for the eigenstates obey the eigenvalue equation 

[U,a]^uja (115) 

For the longitudinal modes, we can take 



a — I d X e 



'3^ ^~ik-x 



cill + c^a + / (iil(c3 J-|_ + CiJ-) 
Using this in the eigenvalue equation, we find 

^i^ f dn Q ■ eS'^Hcs - Ci) 
i k ■ Q 



UJ Ci 



C3 



|fc| {uj + k-Q) 



C\ 



(116) 

(117a) 
(1176) 



and C2 = 0, C4 = c-i{-Q). Using Eq.(117b) in Eq.(117a), we find that a nonzero 
solution requires 

ik-Qf 



UJ / rffJ 

27r./ |/c|2(a;-/c.Q) 







(118) 



This gives the (w, fc)-relation for the longitudinal modes. The creation operator is 
then 



aL(/c) = n(fc)--^ \d^ 
\kV 



k-Q 



-J+{k) 



k-Q 



-J-ik) 



(119) 



(uj + k-Q) (uj-k-Q) 

A similar analysis can be done for the transverse modes, the creation operator being 



arik) = n^ik) + V-^ax{k) - i I dn 



^'^^ -J+{k)- ^'t^\ j-(k) 



(uj + k-Q) 



[LO-k-Q) 



(120) 



So far we have considered plasma waves in Abelian subgroups. One can also 
construct more general solutions ^^. Consider fields of the form A^^ = AJ^(p • x), i.e., 
the field depends only on the combination s ^ p-x where p^ is a timelike vector. The 
zero-curvature condition (75) can be solved as 



a+ = ^A-Q 



I p-Q 
P-Q' 



where A = 44- Writing A^ = e^(p)/ia(s)(-zt"), we find the equations of motion. 



ds 



P^P^^h +nlh3 = o 

p'^Pf^hi + fl'^hi + {hi + hl)hi = 



(121) 



(122a) 
(1225) 



(We have taken the gauge group to be SU{2) for simplicity; there are also equations 
with cyclic permutations of the labels.) Define 



fc^ = 



P^f7 



k^ =^n^ 



(123) 



for transverse and longitudinal cases. The longitudinal mode has a solution of the 
form e**"'^; the relation between fc" and k is given by k^^k^ = 0|, which is the same 
relation as for the Abelian waves. For the transverse modes, the (fc",fc)-relation is 
again the same as for the Abelian case; the amplitudes depend on r = A; • x and obey 
the nonlinear equations 

/a + [1 + {eabfbfifa = 0, a, 6 - 1, 2 (124) 

where now / = ^. These have been further studied in ref.[29]. 

11. Magnetic screening and magnetic mass 

Hard thermal loops and the effective action are important in setting up the 
resummed perturbation theory. Now, r[A] is also an electric mass term for gluons 
and properly incorporating r[A] in any calculation eliminates some of the infrared 
singularities. However, there would still remain some singularities since the static 
magnetic interactions are not screened. It is possible to introduce an infrared cutoff 
by hand to screen the magnetic interactions and make perturbative calculations 
well defined. Interestingly, a gauge invariant magnetic screening term can also be 
constructed using Chern-Simons related techniques ^°. Of course, in principle, one 
should not have to introduce such a cutoff by hand. It is generally believed that for 
QCD at high temperatures there is a magnetic mass term which screens the static 
magnetic interactions (or more generally magnetic fields with spacelike momenta). 
One way to understand how this might happen is as follows ^^ In the standard 
imaginary-time formalism for equilibrium statistical mechanics, bosonic fields are 
periodic in the imaginary time t with period 1/T, i.e., (f){x,T + 1/T) = </>(:?, t). The 
propagator has the form 

ujn = 27rnT, n = 0, ±1,±2,..., are the Matsubara frequencies. At high temperatures, 
because of the structure of the denominator, we expect only T = mode to be 
important. There is no r-dependence for this mode, and thus we expect QCD to 
reduce to three-dimensional QCD. The coupling constant of this reduced theory 
is \J g^T. For QCD in three dimensions we expect a mass gap (~ g^T) and this is 
effectively the magnetic mass of the high temperature four-dimensional QCD. The 
magnetic mass term is expected to be gauge invariant and parity even, although not 
necessarily local. It must be Lorentz invariant if we include the overall motion of the 
plasma. Also we expect such a term to be relevant only for the spatial components 
of the gauge potential. A term which has all these properties is given by 

f [A] = -M2 /" d^ K[An, An] (126) 

where An = ^Aim, An = ^A^fii and 

Ui = {— cos cos (fi — i simp, — cos simp + i cos (p, sin9 ) (127) 



(This is in the rest frame of the plasma.) Notice that n^ is a complex three- 
dimensional null vector; it takes over the role of the null vector Q^ of the electric 
mass term, f involves, in the rest frame, only the spatial components of A^, as 
expected for a magnetic mass term. /(A„), I{An) are defined using n ■ x and n ■ x 
in place of z,z in Eq.(64). Thus K in Eq.(126) is K[^A ■Q,^A- Q'] of Eq.(68) with 

Consider the simplification of the quadratic term in Eq.(126). Using 



dD, riiTij = — 5ij 



3 

k ■ h Stt 

ail n^Uj = — - 

k ■ n 3 



ki Kj I X * J 



(128) 



we find 



Thus f does give screening of transverse magnetic interactions, with a screening 
mass M. 

The higher order terms in f can also be evaluated in a straightforward fashion, 
noting that the basic change is replacing Q^ by {0,ni) and Q'^ by (0,fii). 

Of course, there can be other ways of representing magnetic screening. What is 
notable about the above way of writing the magnetic screening effects is its evident 
kinship to the electric mass term and the Chern-Simons eikonal. Actually, in an 
arbitrary frame, both r and f look the same, except for the choice of different 
orbits of the Lorentz group. Of course, we do not know what M is; no calculation 
even with summations of an infinity of diagrams exists at this point. There are 
some estimates based on lattice and other approaches to three-dimensional QCD. 
Even without a value for M, one can use f as a gauge-invariant infrared cutoff in 
perturbative caluclations. One can also try to estimate M by some self-consistent 
method. Recall that for the Nambu-Jona Lasinio model ^^, once we know the spinor 
structure of the mass term, viz., mqq, we can write 

C — Cnjl + rnqq — mqq = Cq — mqq (130) 

We then calculate the one-loop effective action using £o, taking the subtracted 
term -mqq formally as a one-loop term. The vanishing of the correction to the 
mass then becomes the gap equation determining m. A similar calculation can 
be attempted here. The trouble is that one cannot really stop with the one-loop 
term; contributions from all loops are generally significant. It is possible that the 
self-consistency equations will simplify in some cases, such as perhaps the planar 
diagram approximation. 
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